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Abstract. We investigate the interplay of onsite Coulomb repulsion and various 
mechanisms breaking the fourfold lattice symmetry in a three-band model for 
the iron planes of iron-based superconductors. Using cluster-perturbation theory 
allows us to locally break the symmetry between the x- and {/-directions without 
imposing long-range magnetic order. Previously investigated anisotropic magnetic 
couplings are compared to an orbital-ordering field and anisotropic hoppings. We 
find that all three mechanisms for a broken rotational symmetry lead to similar 
signatures once onsite interactions are strong enough to bring the system close to 
a spin-density wave. The band distortions near the Fermi level are independent 
of differences between the total densities found in xz and yz orbitals. 
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1. Introduction 

Like cuprate superconductors, many of their iron-based cousins have an antiferromag- 
netic (AFM) phase in the phase diagram near the superconducting phase [TJ [2] . As 
phonons are moreover not believed to be strong enough to explain the relatively high 
transition temperatures of pnictides [3], the AFM interactions have attracted large 
interest [I] . Apart from the metallic character of the spin-density-wave (SDW) phase 
in pnictides - as opposed to the Mott insulator in undoped cuprates - a second differ- 
ence is that the SDW in pnictides breaks the four-fold symmetry of the iron-arsenic 
planes down to a two-fold symmetry. 

This breaking of the four-fold lattice symmetry is seen in the conductivity [5] [6] 
and in angle-resolved photo-emission spectroscopy (ARPES) |S1 [HI HH] even at 
temperatures above the onset of magnetic long-range order. While there is a structural 
phase transition at slightly higher temperature and while the in-plane lattice constants 
thus break the rotational lattice symmetry [llj . the effects in experiments appear 
stronger than can be explained by slightly different lattice constants. Additional 
symmetry-breaking of the electronic degrees of freedom has thus been suggested to 
be involved, especially (i) a breaking of the orbital symmetry, i.e., a lifting of the 
degeneracy of the xz and yz orbitals [T2j Q2J [14], and (ii) a nematic phase of the 
spin degree of freedom [T5l [16] . In the latter scenario, magnetic correlations already 
select the preferred ordering vector, e.g., (ir, 0) over the equivalent (0,7r), but do 
not yet establish long-range magnetic order (XTl [HI HU [19] . A related picture involves 
clusters with short-range magnetic order whose AFM direction is pinned by the lattice 
anisotropy [20] . 
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One problem in deciding between the scenarios is precisely that they all break 
the same symmetry, which implies that as soon as one of the breaks rotational 
symmetry, the symmetry in the other degrees of freedom will be broken as well. 
Lattice [5T] [55] [T5] and orbital [T5] US] degrees of freedom thus strongly interact 
with the spin. Nevertheless, identifying the signatures of various modes of symmetry 
breaking may help to elucidate the most important effects. We calculate here the 
spectral density of the three-orbital model, see Sec. El where the four-fold rotational 
lattice symmetry is broken by (i) orbital order, (ii) anisotropic hoppings and (iii) 
short-range AFM interactions. In order to be able to keep the magnetic interactions 
short-range and in order to include onsite Coulomb interactions, we employ cluster- 
perturbation theory, extending an earlier study [23] of anisotropy driven by short- 
range AFM and ferromagnetic (FM) interactions. Anisotropic magnetic interactions 
can distort the spectral density in agreement with ARPES both in otherwise 
noninteracting three- and four-band models and for a regime close to the SDW 
phase [53], orbital order in the non- interacting model was discussed in [TJ]. Here, 
we compare the two scenarios near the SDW in more detail and also compare them to 
symmetry breaking via anisotropic hoppings. While distortions induced by the latter 
cannot even qualitatively be reconciled with ARPES for the noninteracting model, we 
are going to see that results become more realistic near the SDW transition. 

2. Model and Method 

As we want to use exact diagonalization to solve a fully interacting model on a small 
four-site cluster, we have to restrict the model to at most three bands. We use a 
variant [53] of the three-band model proposed in [5J] that gives a better fit of the 
Fermi surface. The unit cell of the Fe-As plane contains two iron and two arsenic ions, 
however, an internal symmetry of this unit cell allows us to write the Hamiltonian in 
terms of a one-iron unit cell (25] [26] as long as we consider isolated planes, as is done 
here. In momentum space, the tight-binding Hamiltonian can then be written in terms 
of a pseudo-crystal momentum k, which is defined as k = k for the xz jyz orbitals and 
k = k + (-7T, 7r) for the xy orbital and which is taken from the enlarged Brillouin zone 
corresponding to a one-iron unit cell. When translating the spectral density A(k, w) 
from pseudo-crystal momentum k to 'lab-space' momentum k, spectral weight of xy 
character is shifted by (tt,tt) with respect to xz and yz weight [571 [HI [53] , 
The k-dependent Hamiltonian in orbital space can then be written as 

H TB (k)= V T^(k)4 d- , (1) 

where d- (dj- ) annihilates (creates) an electron with pseudo-crystal momentum 

k and spin a in orbital v. The orbital indices u, (j, = 1, 2, 3 refer to the xz, yz and xy 
states of the iron 3d manifold, respectively. The T^^(k) = T^' v (k x ,k y ) defining the 
hoppings are given by 

T n/22 = 2t 2/1 cos k x + 2t 1/2 cos k y 

+ At 3 cos k x cos k y (2) 
± 2t 11 (cos2fc x — cos2k y ) 
+ 4ti2 cos 2k x cos 2k y , 
T 33 = A xy + 2t 5 (cosk x + cos k y ) (3) 
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its cos k x cos k y + 2t g (cos 2k x + cos 2k y ) 

4tio(cos 2k x cos k y + cos k x cos 2k y ), 

= Ati sin fc x sin k y , (4) 

= 2it 7 sin fca; + 4zi 8 sin k x cos fc a , (5) 

= 2it 7 sin /c y + Aits sin fc y cos k x , (6) 

where a bar denotes the complex conjugate. Hopping parameters are the same as 
in [23]: h = -0.08, t 2 = 0.1825, i 3 = 0.08375, t 4 = -0.03, t 5 = 0.15, t 6 = 0.15, 
t 7 = -0.12, t 8 = -t 7 /2 = 0.06, iio = -0.024, t n = -0.01, t 12 = 0.0275, = 0.75. 
The chemical potential fi depends on the interaction terms and is chosen to fix the 
filling at 4 electrons per site, for non-interacting electrons, we find fi — 0.47. All 
energies are given in eV. These bands are only an approximation and a three-band 
model may not be detailed enough to capture material-dependent properties. 

We also include the onsite Coulomb interactions including Hund's rule coupling 
and pair hopping [28, 29 . While the couplings for the xz/yz doublet can in principle 
differ from the ones involving the xy orbital, we choose them here to be the same and 
employ the symmetric relations U = U' + 2 J for simplicity. 

Hi n t = U ni, a[t ni, Q 4 + (U' - J/2) n iiQ n i)/3 

i.a i,a</3 

- 2 J S i« ■ S i/3 ( 7 ) 

i,a</3 
i,a</3 

where a, /? denote the orbital and Si, Q («i, Q ) is the spin (electronic density) in orbital a 
at site i. The electron-spin operators are given as usual by Si^ = \ ^2 ss i dl us cr ss ,d il/s ,, 
where er = (a x ,a y ,a z ) is the vector of Pauli matrices. We choose here U = 1.02 eV 
and J = J7/4, because the system is then very close to the SDW transition. The 
interacting system has been found to be very susceptible to magnetic interactions 
breaking rotational symmetry [23j . and we thus concentrate on this regime. 

Similar to the approach chosen in [23) . we explicitly break rotational symmetry 
by Heisenberg couplings that act only locally within the small cluster directly solved 
with exact diagonalization (see method below). Extending the analysis in [53], the 
coupling J x along the a;-direction and J y along y can have different magnitudes and 
the same or different signs: 

-^Heis = ^ x ' + Jy J]] Si, M • Sj :l/ , (8) 

<ij>l|as <i,j>l|y 

where fi, v denote orbitals and (i, j) || x/y nearest-neighbour (NN) bonds along the x 
and y directions. For J x / y > 0, the coupling is AFM. 

We compare this magnetic symmetry breaking addressing the anisotropy of the 
magnetic state to a an orbital symmetry breaking addressing the symmetry between 
the xz and yz orbitals. The latter is implemented as a difference in onsite energy 

^orb = A J2( n i,yz - ni, xz ), (9) 

i 

where A > favours occupation of the xz orbital, as proposed as an explanation of 
the anisotropic spectral density [14] . Finally, as a third way to induce an anisotropy, 
we make hoppings along one lattice direction five to ten percent larger. 
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Figure 1. Spectral density A(k,uj) for onsite interactions near the SDW 
transition (U = 1.02 eV, J = C//4) and a phenomcnological field A = 0.1 eV, 
sec breaking the orbital xz/yz symmetry and favouring the xz orbital. Solid 
lines give the noninteracting bands in terms of the pseudo-crystal momentum 
k; shading gives the spectral weight of the interacting system for 'lab-space' 
momentum k, xy weight is thus shifted by (vr,7r) with respect to k. In the online 
colour figure, red, blue and green shading illustrate spectral weight in the xz, yz, 
and xy orbitals respectively. 

Following [23], this Hamiltonian is treated with the variational cluster approach. 
This method allows to include correlations within a small cluster (4 sites for the 
three-orbital model discussed here) , which is solved almost exactly with Lanczos exact 
diagonalization. Hopping between the clusters is then included as a perturbation [301 
I31j . Long-range order can be treated by optimizing the grand potential with respect 
to fictitious ordering fields [32J [33J , the SDW state of a two-band model for pnictides 
has been studied with this approach [34] [35]. In fact, all parameters of the one-particle 
part of the Hamiltonian can in principle be optimized, we found here that optimizing 
an overall fictitious chemical potential is necessary near the SDW transition to obtain 
a stable solution. 

3. Results 

It has been shown that both orbital ordering |14| and anisotropic magnetic 
correlations [53J can reproduce the band distortions in a manner broadly consistent 
with ARPES, i.e., the yz states at X = (ir,Q) move to higher energies than the xz 
states at Y = (0, tt). Differences between these two ways of breaking rotational lattice 
symmetry mostly affect states around T = (0,0), where an orbital energy splitting 
induces stronger distortions than the magnetic fluctuations in the absence of onsite 
interactions |23| . As can be seen by comparing figures [1] and [2] this remains true when 
onsite interactions push the system closer to the SDW transition. One may also note 
that despite the sizable onsite interactions, one still needs an orbital energy splitting 
A = 0.1 eV to induce the band anisotropy between X and Y as seen in figure Q] The 
anisotropy in the final band structure is here not even extreme, with the yz states 
not quite reaching the Fermi level, and A of a similar order of magnitude induces 
comparable distortions in non- interacting bands [HI d3J . This is in contrast to short- 
range magnetic couplings that become more effective at distorting the bands when 
onsite interactions are switched on |23j . Onsite interactions thus enhance the tendency 
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Figure 2. As figure [T] but without orbital energy splitting. Instead, there are 
anisotropic AFM couplings, sec (BJ, acting within the directly solved four-site 
cluster, i.e., on a very short distance only. J x = 0.04 eV 3> Jy = 0.01 eV 




Figure 3. As figure [T] but without orbital energy splitting. Instead, there is 
a phcnomcnological difference in hopping parameter t^: it is 10 % larger along 
y-direction. 

to short-range magnetic correlations but do not appear to strengthen tendencies to a 
direct orbital energy splitting in a comparable manner. We did not find spontaneous 
symmetry breaking between the xz and yz orbitals, which can be investigated in the 
variational cluster-perturbation theory |321 133) . 

Figure [3] illustrates an alternative way to induce an anisotropy: The hopping 
parameter t% is chosen 10 % larger along the y-direction, a rather larger hopping 
anisotropy For no n- interacting electrons (not shown), this simply increases the 
dispersion somewhat, but does not raise the states at Y. As can be seen in figure [3J 
the interplay of the hopping anisotropy with onsite interactions, which favour a (it, 0) 
or (0, 7r) SDW, distorts the bands: The yz band going from r to X has become very 
incoherent, and the only remaining coherent states at X are above the Fermi level. 
It should be noted that the effect of an anisotropy on various hopping parameters is 
not consistent: It is largest for t2, which is the larger NN hopping entering the kinetic 
energy of the xz and yz orbitals, see (J5J). 10 % anisotropy of the other NN hoppings 
ti, ts and tf leads to far smaller effects (not shown). 

Finally, we take a closer look at the short-range magnetic interactions that were 
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shown to induce band anisotropies in [53] , where the strength of the AFM (along x) 
and FM (along y) couplings were chosen with opposite sign, but of equal strength. By 
varying the two parameters independently, we found that it is not necessary to have 
one FM and one AFM direction. For example, figure [5] shows a spectrum obtained 
for a case where both couplings are AFM, but the one along x is much stronger 
(J x = 0.04 cV) than the one along y (J y = 0.01 cV). The AFM couplings act only 
within the directly solved cluster, i.e., they favour AFM bonds more along x than along 
y. The clusters are coupled within cluster-perturbation theory only via the kinetic 
energy, i.e., there is no long-range magnetic order. The anisotropic signatures seen in 
the spectral density are comparable to those found for J x = —J y = 0.015 eV [23]. As 
the order parameter for nematic order is of higher order, its spontaneous symmetry 
breaking can not be studied in cluster-perturbation theory. 

4. Discussion and Conclusions 

We investigated how various mechanisms of breaking the four-fold lattice symmetry 
of a three-orbital model for Fe-As planes manifest themselves in the spectral density. 
Onsite interactions bring the system here close to a SDW transition, where short- 
range magnetic couplings have previously been found to be more effective at breaking 
rotational symmetry than in the non-interacting system. We find that an orbital 
energy splitting and anisotropic hoppings can lead to qualitatively similar features in 
the interacting bands as short-range magnetic correlations: the states at X can be 
found at higher energies than those at Y and can even move up to the Fermi level, 
see also previous studies [UJOSI]. The band/orbital anisotropy near the Fermi surface 
can be similarly pronounced both in cases with a strong total orbital polarization 
n IZ — n yz « 0.38 (found for A = 0.1 eV as in figure [l} and with nearly vanishing 
polarization n xz — n yz m 0.02 (found for the anisotropic AFM couplings as in figure [2J. 
A splitting of w 100 meV between the xz at Y and the yz states at X, which moves 
the latter close to or just above the chemical potential, can be induced by (i) an orbital 
energy splitting of A = 0.1 eV, (ii) 10 % anisotropy in the hopping parameter ti and 
(iii) a magnetic anisotropy of 30 meV between the x and y directions. 

There are differences in the results obtained in the three scenarios: In the case of 
an orbital splitting, the changes around T = (0, 0) are more pronounced than in cither 
the magnetic scenario of the scenario with anisotropic hopping. In the latter case, the 
spectra appear less coherent even near the Fermi level, and the yz band going from 
r to X almost disappears except for states very close to X that are above the Fermi 
level. Comparing to ARPES [7] |8j E] [10] , this lack of coherence does not appear to be 
in good agreement. Distinction between the other two scenarios is more difficult and 
has previously been discussed for an orbital energy difference in the non-interacting 
model |23j . The conclusions drawn there remain valid for the interacting model: 
ARPES data look more consistent with slighter changes around T than those arising 
from an orbital energy difference large enough to raise the states at X to the Fermi 
level. However, it has to be stressed that both our model (including only the three most 
important orbitals) and our method (where the impact of onsite correlations is only 
treated exactly within a very small four-site cluster) imply substantial approximations. 
The most robust conclusion to be drawn might thus be that a decision between the 
scenarios based on experimental ARPES data remains difficult. Finally, it is also 
thinkable that different mechanisms are driving the breaking of the fourfold lattice 
symmetry in various compounds of the pnictide family. 
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